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Abstract 

We prove that the empirical spectral distribution of a (di,dfl)-biregular, bipartite 
random graph, under certain conditions, converges to a symmetrization of the Marcenko- 
Pastur distribution of random matrix theory. This convergence is not only global (on 
fixed- length intervals) but also local (on intervals of increasingly smaller length). Our 
method parallels the one used previously by Dumitriu and Pal (2012). 



> 

o 

o\ . 

tJ- . 1 Motivation 

(^ . In classical random matrix theory there are two basic types of symmetric ensembles: Wigner 

^^ I matrices and Wishart-like ones. There is a simple parallel between them, roughly expressed 

in the following way. Given a non-symmetric real matrix G, there are two natural ways to 

construct from it a symmetric matrix: if G is square, one way is to consider its symmetric part 

^ I A = 2 ' ^^^ other way works for rectangular matrices, too, and consists of multiplying it 

H I by its transpose: W = GG . If one starts with a random matrix G with i.i.d. entries of norm 

and variance 1, if G is square, the first symmetrization yields a Wigner matrix; the second 

symmetrization yields a Wishart-like matrix (it is Wishart, more precisely central Wishart, 

if G consists of standard normal variables; we call it Wishart-like otherwise). 

The spectra and eigenvectors of Wigner and Wishart-like matrices have been shown to 
exhibit universal behavior: many of their eigenstatistics have limiting distributions which 
are independent of the entry distribution, modulo certain technical conditions (in addition to 
being mean 0, variance 1). These results, along with successive weakenings of the technical 
conditions, are the subject of a recent set of breakthrough papers [TVllt ITVlOl lESYOObf 
lESY09a[ lEPR+lol lERSYlOl lERS+lOJ . 

A natural question for the discrete probability community is whether this universal be- 
havior extends to adjacency matrices of random graphs; we are specifically interested in the 
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case of random regular or semi-regular graphs. Such graphs are known to have very inter- 
esting properties: they are good expanders, some classes are quasi-Ramanujan, they have 
wide spectral gaps and as such they mix rapidly, and are of interest in computer science and 
engineering and in coding theory. 

Ordinary random regular graphs (or d-regular graphs), where every vertex has the same 
degree (which grows as a function of the number of vertices), have been recently investigated 
in |DP12[ IDJPP121 IJohl2l [JPT2l ITVW13[ IBL13J : we aim to extend some of the results to 
bipartite, biregular random graphs, where the two sets of vertices have the property that all 
vertices in the same set have the same degree (also growing with the total number of vertices). 

The question of whether the spectra of random d-regular graphs have the same behavior 
as the spectra of Wigner matrices is non-trivial in nature, since the adjacency matrices of 
regular graphs have strong dependencies, namely, all rows and columns add to the same 
number — their common degree. As such, they are not Wigner; in fact, for d fixed, McKay 
|McK81) showed that the scaled empirical distribution function (or ESD; defined below) of 
random d-regular graphs on n vertices converges in probability (and almost surely, if the 
random graphs are defined on the same probability space) as n — )• oo to the Kesten-McKay 
distribution, which has density 



otherwise. 

This differs from the Wigner matrix case, where the scaled ESD converges in probability to 
the semicircle law, which has density 



otherwise. 

When d is allowed to grow with n, the scaled empirical spectral distribution of the random 
d-regular graph does converge in probability to the semicircle law (see [DP 12] and |TVW13] ) . 
Thus, even if for d fixed they are rather different, the spectra of d-regular graphs with d and 
n growing to infinity are similar to the spectra of large Wigner matrices. 

Motivated by a question asked by Babak Hassibi, we study here the spectra of bipartite, 
biregular random graphs. We find that their spectra have similar behavior to the spectra 
of Wishart-like matrices; at first glance, this may appear suprising, but further examination 
reveals linear algebraic reasons why this should be so. Our notation and main results are pre- 
sented in Section [2j In Sections [3] and HI we prove global and local convergence, respectively, 
of the ESD to its limiting measure. The appendix contains the proofs of some statements on 
the distribution of cycles in biregular bipartite graphs, used in Sections [3] and [H to show that 
our graphs are locally well approximated by trees. 

2 Preliminaries and statements of results 

We assume that graphs do not have loops or parallel edges. The adjacency matrix of a graph 

G is defined as 

1, if z ~ j, 



' U, otherwise. 
Note that A is symmetric, and therefore all of its eigenvalues are real. 



Bipartite graphs are graphs composed of two sets Vi and V2 of vertices, with edges only 
between vertices from Vi and vertices from V2; with proper labehng, their adjacency matrices 
have the special form 

X 
X^ 



B 



where the matrix X is defined by the edges between the two classes of vertices. Note that if 
M and N have sizes m and n, respectively, with m < n, then X is an m x n matrix of Os 
and Is. It is a simple linear algebra result that the non-zero eigenvalues of B come in pairs 
(—A, A), with A > an eigenvalue of XX'^ , and that B has (at least) n — m eigenvalues equal 
toO. 

If G is a bipartite graph with vertex classes L and R, then we say it is (dL,dR)-biregular 
if all the vertices in L have degree d^, and all the vertices in R have degree d/j. For simplicity, 
we will always assume that dn > di (and therefore that \L\ > \R\). 

By a random ((i/,,(iij)-biregular bipartite graph with (m + n) vertices we mean a graph 
selected uniformly from the space of all (d/,, dij)-biregular bipartite graphs with \L\ = m and 
\R\ =n. 

We define the empirical spectral distribution or ESD of a symmetric n x n matrix M to 
be the probability measure ^in on the real numbers given by 



^E*^.. 



fJ-n 

n . 

1=1 

where 5x is the point mass at x and Ai, . . . , A„ are the eigenvalues of M. Note that if M is 
random, // is a random probability measure. 

We say that a sequence ^i,//2, • • • of random probability measures on the real numbers 
converges almost surely to a deterministic probability measure // if as n — >• 00, 

f d^Xn -^ fdfi a.s. 

for all bounded continuous functions /: M — )• M. This is equivalent to the slightly different 
statement that with probability one, fin converges weakly to fi. 

Following combinatorialists' conventions, we will often refer to a random graph G when we 
really mean a sequence of random graphs with an increasing number of vertices (depending 
on m and n). We assume that all of these random graphs are defined on the same probability 
space, but we make no assumptions about their joint distribution; all of our results hold for 
any arbitrary joint distribution, so long as the marginal distributions are as described. Many 
variables that we mention implicitly depend on n, and asymptotic expressions O(-) or o(-) 
reflect behavior as ?tt,, n — )• 00. We will occasionally use the notation OAi') to indicate that 
the constant in the big-0 expression depends on some other constant A. 

It is a standard result in random matrix theory that if X is an tti x n random matrix whose 
entries are iid real random variables with mean zero and variance one, and m/n converges to 
a finite limit, then the ESD of -X'^X converges to the Marcenko-Pastur law (see [BSIOJ . for 
example). We show here an analogous result: 

Theorem 1. Let G be a random {dL,dR)-biregular bipartite graph on m + n vertices, with 



the following conditions on the growth of di and dji: 

lim dfi = oo, (1) 

for any fixed e > 0, dji = o{n^), (2) 

Let ^ = ( ^T ) be the adjacency matrix of G (under proper labeling). Then as n ^f oo, the 
ESD of —X'^X converges almost surely to the Marcenko-Pastur law with ratio y^^ ■ This 
distribution is supported on [a^,6^] and is given on that interval by the density 

y 



P{^) = 2 — \/(^^ - x){x -a^), 

where a = 1 — y~^''^ and 6 = 1 + y~^''^ . 

Theorem [1] agrees with the results of |MS03] . in which Mizuno and Sato derive the hmiting 
distribution of the eigenvalues of a sequence of deterministic biregular graphs with girths 
growing to infinity. They do so using the Ihara zeta function, and they express their result 
in a different form from ours, which turns out to be equivalent. It should be noted that our 
result is much stronger than theirs; as we will see, most bipartite biregular graphs have small 
girth, but this does not affect the convergence of the ESD. 

As mentioned above, the ESD of d^ XX is the distribution of the squares of the non- 

— 1/2 

trivial eigenvalues of dj^ A. We can thus find the limiting distribution for the ESD of this 
matrix as well: 

— 1/2 

Corollary 2. The ESD of dj^ A converges almost surely to the distribution p, supported on 
[—6, —a] U [a, b] and given on that set by the density 

2 , 2m , y 



p{x^)\x\= \ J{b^-x^){x^-a^l (4) 

along with a point mass of ^-^ at 0. 

It is also known that when d — >• oo, the ESD of random d-regular graphs converges to 
the semicircle law on short scales (see [DP 12) , |TVW13 j ) . In Section [H we prove this for the 
biregular case, under slightly different conditions on the growth of dji: 

Theorem 3. Let G be a random (di, dR)-biregular bipartite graph on m+n vertices satisfying 
([TJl-Q, as well as the more stringent condition dfi = exp(o(l)-y/logn). Fix e > 0. Let A be 
the adjacency matrix of G and jin be the ESD of [dji — l)~^'^A, and let jjl be the limiting ESD 
of Corollary l^ There exists a constant C^ such that for all sufficiently large n and 5 > 0, for 
any interval I C M avoiding [— e,e] and with length \I\ > max {2r], r] / {—5 log 6)) , it holds that 

|/i„(/)-/i(/)| <6C,\I\ 

with probability 1 — o(l/n). The quantity r], which gives the minimum length of an interval L 
that we consider, is given by the following series of definitions: 

. / logn 
\9{log dRY 
r = e^/% 

^^l/2_^-l/2_ 



This theorem is obscured by the technicahties in its statement, so we give some discussion 
of its meaning. Corollary [2] only gives information on ^n{I) for fixed |/|. Theorem [3l on the 
other hand, allows |/| to shrink as r/ does. We note that r/ ~ 1/a. 

We restricted our intervals / away from to avoid complications caused by the point mass 
that /i has when dpi/di — )• y > 1. Since the support of ^ except for this mass is bounded 
away from 0, this restriction costs us nothing. 

To prove our results, we use the moment method along with Stieltjes transforms, combined 
with a careful examination of the local structure of the bipartite, biregular graph. Along 
the way, we need to adapt some of the results proved by McKay, Wormald, and Wysocka 
|MWW04] for random d-regular graphs to our {dL,dfi) biregular random graphs. We give 
these results below. The method used follows |MWWd4] very closely, which is why we relegate 
the proofs to the appendix. 

Let G be a random (d^, (i/j)-biregular bipartite graph on ?n, + n vertices. Assume d^ < dji, 
and let a = dji/dL. As always, all of these variables depend on n, and any expressions O(-) 
or o(-) reflect behavior as n — )■ oo. We assume that a converges to a finite value as n — t- oo 
(this assumption is also necessary in the case of Wishart matrices). Here and throughout 
the paper, "cycle" always refers to a simple cycle, with no repeated vertices. Let Xr denote 
the number of cycles of length 2r in G. (Note that as G is bipartite, its cycles all have even 
length.) 

Proposition 4. Let 

{dL - lYidR - ly 



fi. 



r 



2r 
If dji = o{n), r = O(logn), and rdji = o{n), then 



Var[X^] =Hr[l + 



dj^jra^'-^ + g-'-dR) 



n 



In Section [3l we use this proposition to show that with high probability, our biregular 
bipartite graph G is locally well approximated by a tree. This allows us to approximate the 
traces of the adjacency matrix of G, thus computing the moments of its ESD. Finally, we 
refine these results in Section U] to prove local convergence on vanishing-length intervals. To 
prove this theorem, we give estimates on the rate of convergence of the Stieltjes transform of 
the ESD, the same approach used in [DP 12] and [TVllj . 

3 Global convergence to the Marcenko-Pastur law 

To find the limiting ESD of a biregular bipartite graph G, we will first show that in a sense 
that we will make precise, most neighborhoods in G have no cycles and are trees. This will 
allow us to estimate the traces of the adjacency matrix of G, and we will find the limit of 
these as G grows with a combinatorial argument. 

For this entire section, let G be a random (d^, (i/j)-biregular bipartite graph on n + m 
vertices, and assume that conditions ([I])-® on the growth of dj^ and dji hold. As before, let 
a = -^. Condition ([2]) is more restrictive than the condition dji = o{n) of Proposition Hj it 
implies that for any fixed r and e > 0, we have d^ = o{n''). 

We make precise the property of G being locally a tree in the following lemma: 



Lemma 5. Let r he fixed, and let r he the set of vertices in G whose r -neighborhoods contain 
no cycles. Then, if da satisfies ([I]) and ([2]), 



1 



> n 



-1/4 



>(n-5/4) 



n + m 

This is the same statement proven in |DP12| for regular graphs, and using Proposition [J] 
we can prove it in the same way. 

Proof. If a vertex is not in r, then for some s < r there exists a 2s-cycle within r — s of 
the vertex. The size of all (r — s)-neighborhoods of 2s-cycles hence serves as a bound on 
the number of "bad" vertices. For any given 2s-cycle, the size of its (r — s)-neighborhood is 
bounded by 2s{dfi — l)*"^*. If we define 



A^: 



Y,MdR-ir-''Xs, 



s=2 



then this gives us the bound n + m —\t\ < N*. 

Now we compute E[A'^] and Var[A^*]. Using our expression for E[Xr] from Proposition [H 

r 

E[iv;] = Y, MdR - ly-'idL - lyidR - lyoii) 

s=2 

r 

= idR-iYY.oiidL-ir) 

s=2 

= o{{dL-indR-iY) = o{dj^) 

To compute the variance, first notice that ([2]) implies that Var[Xs] = /is(l + o(l)). By 
C auchy- S chwar z , 



Var[Af;] < r^4s2(d^ _ lf'-^'Yar[Xs] 

s=2 

r 

<r^4s2(d^-l)2-2s^, (1 + 0(1)) 



s=2 



idn - ir E ' 



{dL - lYidR - ly 



s=2 



(dR - 1)2« 

r 

<ridR-iril+oil))Y,s 
The rest of the lemma follows from Markov's inequality: 



(1 + 0(1)) 



T 



n + m 



> n 



-1/4 



P [n + m - |r| > (1 + a)n^/^] 



<P[iV; > (1 + a) 
Var[iV;] + E[iV;] 



3/41 



*12 



< 



(l + a)2^3/2 



O 



3/2 



n 



o[n 



-5/4^ 



D 



This result shows that there are few "bad" vertices. It easily follows that this is true 
within the left and the right vertex classes of G as well. 

Corollary 6. Let tl and th he the number of vertices in the left and right classes of G, 
respectively, with acyclic r -neighborhoods. Then 



P 
P 



m - \tl\ 
n + m 

' n- \tr\ 
n + m 






o(n"5/4). 



Proof Note that 1 

»-|Tfl| 



\r\ 
n+rn 



n-\-m 



> C. 



^^J^ + ^^^, SO 1 - ;^ > c whenever ^^^ > c or 

D 



Let I3k{r,a'^) be the A;th moment of the Marcenko-Pastur law with ratio r and scaling 
factor a'^ as defined in [BSlOj . 

— 1/2 

Proposition 7. Let A be the adjacency matrix of G, and let fin be the ESD of dp^ ' A. 
Recalling that y = lim„_j.oo «; 



^2k+l 



dfin{x) — )• a.s., 



X^^d^XrAx) 



l + a 



Pk{y \l) a.s. 



as n —)■ oo. 



Proof. Consider the infinite (di, (i/j)-biregular tree. Let Br denote the number of closed walks 
of length r on this tree, starting from some fixed vertex of degree di, and let Cr denote the 
number of closed walks of length r starting from a vertex of degree dji. Note that as dL and 
dR depend on n, so do Br and Gr- 

First, we formulate the rth moment of //„ in terms of Br and Gr- 



x"" dflriix) 



d 



-r/2 



R 



n + m 



E ^''(^- 



v£V{G) 



The quantity A^{v, v) is the number of closed walks of length r from v in G. For v G tl, this 
is Br, and for v S tr, this is Gr. For v ^ t, we can use the bound A'^(v,v) < d^. Hence we 
can bound the rth moment of /x„ by 



d 



-r/2 



R 



n + m 



{\TL\Br + \TR\Gr) < / x"" dHn{x) 



< 



d 



-r/2 



R 



n + m 



[mBr + nGr + {n + m — |T|)(i^) . 



Define 



Q-n — (^R 



r/2 



m 



n-y^)Br+^ 

n + m / \n + m 



n 



n 



-V^V. 



^R 



r/2 I rnBr + nGr 



n + m 



+ 



n 



-1/4 



n + m 



-d\ 



R 



Now 

an < /"x'-dMn(x)l > P \l^—M < n-V4 ^nd ^^^M < n-V4 
J J \_ n + m n + m 

= 1 - o(n-5/^), 

and in the same way, P \^^ x^ djinix) < 6n] > 1 — o(n~^'^). By the Borel-Cantehi lemma, it 
holds almost surely that an < J x'^ dfin{x) < bn for all but finitely many n. If we show that 
an and bn converge to a common limit (which we will do next), it will follow that f x^ d^n{x) 
converges to this limit almost surely. 

To find the limits of a^ and 5„ as n ^- oo, we first note that n~^'^(i^ — t- by ([2]). Since 
Br,Cr < d^, this also implies that n~^'^Br — )• and n~^'^Cr — )• 0. So, it suffices to show 
that 



d 



-(2k+l)/2 



R 



n + m 



imB2k+i + nC2k+i) -^ 0, (5) 



-^{mB2k + nC2k) -^ -^h{y-\ 1). (6) 

n + m \ + a 

Equation ([5]) is trivial, since B2k+i = C2k+i = 0. To prove ^, we introduce the Narayana 
numbers (see |Sta991 p. 237]), defined as 

N{k,a) 



a + 1 \aj \ a 

The moments of the Marcenko-Pastur distribution can be given in terms of these numbers 
[BSlOl Lemma 3.1]: 

fc-i 
/3k{y-\l) = Y.y-''N{k,r). (7) 

r=0 

We will give a combinatorial argument to relate the closed walks on the tree to the Narayana 
numbers. We mention that another approach to proving ([6]) is to calculate Bj. and Cr using 
the spectral density of the infinite (dL,d/j)-biregular tree, as calculated in |GM88t (5.7)]. 

A Motzkin path of length 2k is a lattice path that starts at (0,0), ends at (2/c,0), and 
stays above the x-axis; each step can be a rise (/^), a fall (\), or a level step (—)•). An 
alternating Motzkin path is a Motzkin path that rises only on even steps and that falls only 
on odd steps. See Figured] for an example of the five alternating Motzkin paths of length 6. 

The alternating Motzkin paths have the following connection to the Narayana numbers: 

Lemma 8 (Lemma 6.1.7, |Dum03j ) . The number of alternating Motzkin paths of length 2k 
with exactly a rises is N(k,a). 

We relate the Narayana numbers to the walks on a tree by the following two lemmas. 
A ballot sequence of length 2k is a sequence xi, . . . ,X2k of I's and — I's such that all partial 
sums xi + ■ ■ ■ + Xj are nonnegative. 

Lemma 9. The number of ballot sequences of length 2k with a I's at even locations and k — a 
I's at odd locations is N{k,a). 




. _L 1 . 



Figure 1: The alternating Motzkin paths of length 6. 
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Figure 2: An alternating Motzkin path and its corresponding ballot sequence. 

Proof. We give a bijection between alternating Motzkin paths and ballot sequences. Encode 
the alternating Motzkin path as pi, . . . ,P2A:> where pj is 1, 0, or —1 depending on whether the 
ith. step is rising, level, or falling. Define a sequence by Xi = 2pi + (-1)*"^. We will confirm 
that this is a ballot sequence: each Xi is either 1 or —1; for any j, 

j 

xi + ---+xi = 2{pi +---+Pi)+ x](-l)*~^ 

and both of these terms are nonnegative; and xi + • • • + X2k = 2(pi + • • • + P2k) = 0. So, 
xi, . . . , X2k is a ballot sequence. 

To map back from ballot sequences to alternating Motzkin paths, we let pi = (xi — 
(— 1)*~^)/2. For any even j, 

Pi-\ \- Pj = -^{xi-l h Xj) > 0. 

For any odd j, 



Pi H \-Pj = -^{xi-l \-Xj -1). 

Since xi + ■ ■ ■ + Xj > 1 when j is odd, this expression is also nonnegative. So, our path 
stays above the x-axis. The other properties of being an alternating Motzkin path are easy 
to check. 

This bijection takes alternating Motzkin paths with a rises to ballot sequences with a I's 
at even locations, so the lemma follows from Lemma [HI D 



Lemma 10. 

fc-i 



B2k = Y.idR-irdl~^N{k,a), (8) 

a=0 

fc-1 _ 

C2fc = 5](dL-ir4-'^iV(A:,a), (9) 

a=0 

/or some d^ and (Ir satisfying di — 1 < di < di and dn — 1 < dn < d^. 

Proof. Fix some vertex v in the (d^, d/jj-biregular tree with degree di to serve as the root. 
We will enumerate the closed walks of length 2k starting at v. To any such walk we can 
associate a ballot sequence of length 2k, given by putting a 1 at every step of the walk going 
away from v and a —1 at every step returning toward v. We will count the number of closed 
walks associated to each ballot sequence. 

Fix some ballot sequence, and suppose we are constructing a closed walk associated with 
it. For every 1 in the ballot sequence, our walk must go outward from the root. If the 1 is 
at an even location, we have d/j — 1 choices for where to move; if it is at an odd location, 
we have either cIl oi di — 1, depending on whether we are moving from v or from some 
other vertex. For every —1 in the ballot sequence, our walk must move backward towards 
the root, and there is no choice to be made. So, given a ballot sequence with a rises on even 
steps, the number of closed walks from v associated with that ballot sequence is between 
{dn - l)°'{dL - l)''-"' and {dR - 1)"^^"". Using Lemma M to count the number of ballot 
sequences with r rises on even steps, we obtain ([8]). The same proof starting with a vertex v 
with degree d^ gives us Q. D 

Now we finish the proof of Proposition [7] by computing the limit of 

^ (mB2k + nC2k) 



n + m 



as n — 7- oo. Recall that all of these variables depend on n except for k, which is fixed. By 
Lemma [T0| we can rewrite the above expression as 

_jr^ ^ JdR-iyd-,' ^ r) + ^^ V ^-^^Nik, r). 
n + m^ dl-^ n + m^^ d^-^d\ 

Replacing m/[n + m) and n/{n + m) by a/(l + a) and 1/(1 + a) respectively, and taking the 
limit as n — )• oo yields 

fc-i fc-i 

-^ Y, y''-'N{k, r) + -—Y, y-'^Nik, r). 

Note that this is where we used ([T]). We can simplify this expression to 

2 tX 
—-Y,y-^Nik,r). 

^ r=0 

This is exactly j^:zl3k{y^^, 1) as given in d?]). D 

10 



Proof of TheoremUl Let m„ be the ESD of d^ X^X. As described in the introduction, the 

1/2 1/2 

eigenvalues of dj^ A consist of zbuj for the singular values ui, . . . , (T„ of d^ X, along with 



m — n O's. It follows that 



f k 7 "*- + "-/" 2A- , 
X dun = — / X d^„ 



a + l 



X^'^dlln- 



It follows from Proposition [7] that j x^ dvn ~^ h{y~^ ^ 1) a-s- as n — t- oo. Since the moments 
of Vn converge almost surely to the moments of the Marcenko-Pastur distribution, which is 
supported on a compact interval, Vn converges almost surely to this distribution. D 

To wrap things up, we compute the density of the limiting distribution of //„. 

Proof of Corollary [3 We only need to show that the moments of the measure given by the 
density {2 / (1 + y))p{x'^)\x\ agree with the limits of the moments of //„ found in Proposition [71 
The odd moments of this measure are by symmetry, and the even moments are easily 
computed by integrating and substituting u = x'^. D 

See Figure El for a picture of the limiting distribution for a few values of a. 

4 Convergence on short scales 

We will use the method of [DP12| Section 3] . The basic idea is to use the local approximation 
of our graph as a tree to show that the graph's Stieltjes transform is very close to that of a 
tree. We break it down into the following steps: 

1. Compute the resolvent of the adjacency matrix of a biregular tree, and show its (root, root) 
entries converge quickly to their limits. 

2. Let ?; be a vertex of a deterministic biregular graph G with a large acyclic neighborhood. 
Show that the (u, v) entry of the resolvent of G is close to the corresponding entry of 
the resolvent of a biregular tree, hence bounding the rate of convergence of this entry 
to its limit. 

3. Show that nearly all the vertices of a random biregular graph have a large acyclic 
neighborhood, and use this fact to transfer the estimates of the earlier parts to random 
graphs. 

4. Use the method of [TVIH Lemma 64] to deduce Theorem [3] from the estimate on the 
Stieltjes transform. 

Our task is made slightly more difficult by the need to consider two different (d/,, (i/j)-biregular 
trees: one in which the root has degree di-, and one in which the root has degree dn. 
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Figure 3: This figure depicts the Hmiting density of the ESD of a random (di,dR)-biregular 
bipartite graph with a few different values of a = -^. The spike at denotes a point mass. 
The continuous part of the density is given by ([3]) , and the point mass has size ^^ . Each of 
the left and right spikes are scaled down copies of the Marcenko-Pastur distribution under 
the transformation x i— )• ^/x. When a = 1, the density reduces to that of the semicircle law. 
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4.1 Preliminaries 

The resolvent of a Hermitian matrix A is the matrix R{z) = {A — zl)~^, for '^{z) > 0. The 
resolvent of a graph is the resolvent of its adjacency matrix. The Stieltjes transform of an 
nxn matrix A is the Stieltjes transform of its ESD, the function given on the upper half-plane 
hy siz) = ^tT{A - z)-K 

We will use the following, well known formula for the inverse of a block matrix. 

Proposition 11. Let A and D be nxn matrices of size nxn and nix m, respectively, and 
let B be n X m. Let 

M 
Then 



A-^ + A~^BF~^B^A-^ 
-F-^B^A-^ 



A B 








B^ D 








A-^BF 


-1" 


F = D- 


- B^A~^B 



In this section of the paper, we will define the ratio a hy a = {dn — l)/{dL — 1) rather 
than d^/di. Let Un{z) be the Chebyshev polynomial of the second kind of degree n. We 
define the following shifted Chebyshev polynomial. 

This family of polynomials satisfies the recurrence 

9o(^) = 1 

qn{z) = {z^ - a"^ -\)qn-\[z)-a~^qn-2{z), n > 1, (10) 

which follows by applying the recurrence Un{z) = 2zUn-i{z) — Un~2{z)- 

4.2 Resolvents of trees 

Our aim is to calculate the resolvents of biregular trees. We start, however, by considering 
trees in which each vertex has either d^ — 1 or d^ — 1 children; this means that every vertex 
has degree di or dn except for the root, which has degree d^ — 1 or d/j — 1. 

We define T^iC) to be the tree with depth (^ where the root has d^ — 1 children, its 
children each have dji— 1 children, their children each have di — 1 children, and so on. The 
tree Tl{Q) is a single vertex. We define Tji{Q similarly, but with the root having dji — 1 
children. To determine the adjacency matrices of these trees, we place an ordering on the 
vertices as follows: If C = 0, then there is only one vertex and hence one possible labeling. 
For (" > 0, we will define an ordering inductively. Choose a subtree of the root and list of 
all its vertices in the order already determined for C ~ 1- Then, do this with the remaining 
subtrees of the root. Finally, put the root last. 

Let Hl and Hr be the adjacency matrices of Tl{C) and Tr((^), respectively. We define 

fLiO = {{dR - 1)-'/^Hl - ^)root,root 
^L(C) = ((dK-l)-^/^^L-.Ut,,eaf 

We note that V'l(C) is independent of the particular leaf chosen. We will make use of the 
recursive structure of the trees to calculate these values. 
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Lemma 12. (a) 



qc_{z) +a~^g^_i( 
zqc_{z) 

zq^{z) 



V^l(2C+1) 
<^fl(2C + 1) 



zqt^iz) 






(b) 



^l{20 = i^R{20 
V'l(2C + 1) 



{dR - 1)-^ 



V'ii(2C + 1) 



(d^ - l)-f-V2 



Proof. We will start by showing that ^l{C) ^-iid ^r{C) satisfy the recurrences 

^l(C) = -(^ + «-Vk(C-i))"\ (11) 



^r{C) 



+ ^z.(C-l))"'. 



(12) 



Consider the tree T^ of depth ( and let Hi, ... , Hdj^-i denote the adjacency matrices of the 
subtrees of the root. Using the given ordering for the vertices, we have 



^/d 



R 



1 



Hl{0 - z 



Vd_ 



fcT^^(C-l) 



vdR 



=T^«(C - 1) 



Vd 



fcr^MC-i) 



u 



where u is a column vector representing the children of the root. This vector is {dR — 1)^^'^ 
in the root of each of the subtrees of the root and elsewhere. Using Proposition [11] and 
thinking of the —z in the bottom right corner as a 1 x 1 block, we find 



MO 



dL-l 
dB-1 



miC-r 



which is (jlip . The proof for (|12p is the same. 

Unwinding these recurrences and noting that (fiiO) = ^r{0) = —z~^, we have the fol- 
lowing continued fraction representation of ipiiC)' 



MO 



a 
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Using standard formulas for the evaluation of continued fractions (see |LW08| ) , we find that 



</^l(C) = ^, where 



^2C = 


= Z^2C-l - 


- Mc.- 


-2 


-B2C = 


= ZB2C-I - 


- B2(^- 


-2 


with the initial conditions 








Ao = - 


1, 






Bo = z 


; 





-42C+1 = 2^2C - " ^^2C-1 



-1 



-62^+1 = zB2( — a i?2C-i 



Ai = -z, 

Bi = z^ -a-^. 

We can iterate these recurrences as follows: 

^2C = z{zA2Q-2 - or'^A2C_-'i) - ^2C-2 

= z[zA2(;-2 (-42^-2 + ^2C-4)) " -42^-2 



Z 



2 a-^ - l)^2C-2 - a"^A2^_4. 



Applying this procedure to the vl2(^+i and to the i?2C ^^^ -^2C+i cases give the same result, 
yielding 

A(^ = {z^ - a"-^ - l)Af-i - a'^A(^^2, 
Bc_ = (z^ - a'^ - l)5c-i - ""^-Bc-2. 



It is easily checked using (|T0|) that 

A2C = - (qd^) + a''^qc-i{z)) , ^2c+i = -zgciz), 

B2C = zqciz), B2C+1 = qc+i{z) + qc{z). 

From these expressions and ()12p . it is straightforward to derive the expressions for ipji(2Q 
andv9ij(2C + l). 

To compute V'l(2C) and iPr{2(), we will first show 

MO = -{dR - 1)-^'\l{C)^r{(: - 1), (13) 

V'fl(C) = -{dR - i)-^/Vi?(C)V'L(C - 1). (14) 

Using (jlip . if A is the minor of (d/j — 1)"^"-?//, — z consisting of all but the last row and 
column, we have 

1 



^1,l{C)=[^==Hl-z 

kVO-R-^ / l,root 

-MC){A-'u)i 

-(dR-ir^/^iPLiOi^Ric-i), 
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proving (fT3]) . Equation ([H]) is derived in the same way. By combining these, we get •0l(C) 
{cIr - 1)- VL(C)¥'ii(C - l)MC - 2), whence 



M'iC) = MO){dR - ly'^ n Vl{2j)^r{2j - 1) 

The proof of the expression for Vi?,(2C) is identical, and the expressions for V'lC^C + 1) ^^^ 
tpni^C + 1) follow immediately by applying (fTSl) and (fT4|) . D 

We now turn from these almost regular trees to the real thing. Let Tl{C) be the {dL,dR)- 
biregular tree of depth ( whose root has degree cIl, and let Tr{() be the (dL,dR)-biregular 
tree of depth ( whose root has degree dR. Let H^ and Hr be the adjacency matrices of T^iC) 
and Tr{(), respectively, with vertices ordered as with H^ and Hr. We define 



MO = {{dR - ly'/^HL 



/root, root' 
'root. leaf 



^r{C) = {{dR - 1)-^'^Hr 



^I^r{C) = {{dR 



-1/2 



Hr-z\ 



7 root, root ' 
root, leaf 



Lemma 13. 

^l(2C) = 



V'l(2C) = ^r{K) 



^(9c(^)-d^9C-i(^))' 

{dR - !)-< 



M2C) 



q^{z) + q<^-l{z) 
z{qd^) - di=Tl(~l{^)) 



z{Qd^) - d]^q(-i{z)) 

Proof. Because the root of Hl has di children instead of (i^ — 1 , the methods of Lemma [T2] 
give 



^l(C) 



dL 



dR-l 



^r{C-i) 



Substituting in the value for ipR{C — 1) from Lemma | 
expression for (pr{2Q is derived in the same way. 
The same procedure shows that 



I yields the desired expression. The 



V'l(2C) = -{dR - 1)-i/Vl(2C)Vr(2C - 1), 



D 



and substituting the value from Lemma [12] yields the desired expression. 

We will now bound the rate of convergence of some of these functions to their limits as 
C — >■ oo. First, define the complex function F{z) = z + \/ z^ — 1, with branch cut [0, oo) for 
the square root. Let ^(z) = F{\^Iol{z^ — oT^ — 1)) and r{z) = \w{z)\. We will refer to 'w{z) 



16 



and r(z) as simply w and r. Note that r > 1 for all '^{z) > 0. Using a well known expression 
for the Chebyshev function Un{z) (see [MH03J ). we can expand q({z) as 



qt^{z) = a 



w — w~^ 



(15) 



As we will see, the limits of 93l(C) ^nd (Pr{C) as C — ^ c» are given by the following 
functions, defined on the upper half-plane: 



SL{Z) 



1 g-^/^w-' 

z z 



sr{z) 



1 ai/2^-1 



with branch cut [0,oo) for the square root. We define 

asL{z) + SR{z) a 



s{z) 



1 + a 



1 + a 



z + \ z 



1 — a 



az 



again with branch cut [0, oo) for the square root, and we note that s{z) is the Stieltjes 
transform of the limiting ESD /i of Corollary [2l 

We bound the convergence in terms of r = r{z) and (: 



Lemma 14. 



\ipL{2C) - sl{z)\ < 



\iPRi2C) - sr{z)\ < 



\M2C)\, l^il(2C)l < 



|z|(l-r-2C-2) 
2aV2^-2C 



|z|(l-r-2C-2) 



2a?/2(dfi-i; 



-Cr-C 



|z|(l-r-2C-i) 



Proof. Applying ([15]) to the formula for ipl(2() 

a-i/2 



\ipU20 - sl{z) 



a 



-1/2 



_i w'' — w 


c 


y;C+l — yj 
^-2C-l(l_y;-2) 


C-1 



1 - 'u;-2C-2 



< 



2a-l/2^-2C-l 

|z|(l-r-2C-2)' 



The exact same procedure establishes the corresponding inequality for \ipR(2() — sr{2()\. We 
can similarly compute 



\M2C)\ = l^i?(2C)l 



{dR-l)-<{w-w-^) 



< 



z(it;C+i-u'-C-i) 

2a</^{dR-l)-<r-< 

U|(l-r-2C-2) • 



D 
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4.3 From trees to deterministic graphs 

We now move from trees to graphs with large, acychc neighborhoods. Let G be a deterministic 
{di, dii)-hiregular graph that has a root vertex with an acychc ((" + l)-neighborhood. Let A 
be the adjacency matrix of this graph. Our goal is to show that the resolvent of A is well 
approximated by the resolvent of H^ or Hji. If the root of G has degree di, then we consider 
the error term 

El{0 = {{dn - I)-'/' A - z);^^^ _^ - ^UO, 
and if the root of G has degree d^, then we consider the error term 

EniC) = {{dn - I)-'/' A - z);^^^^^^^ - ^n{C)- 
Lemma 15. We can bound El and En by 

\MC)\'dL{dR-iy^/'\dL-l)^'^/'i 



\El{C)\ < 
\Er{C)\ < 



{dn - l)9(z) 

\MC)\'dRidL-i)^^/'HdR-i)i</^i 

{dR - 1)9(Z) 



Proof. Let Cl and Cr denote the number of vertices of distance C + 1 from the root in the 
biregular trees Tl and Tr, respectively. The proof of [DP121 Lemma 10] goes through gives 
the bounds 



(dR - 1)9(Z) ' 



|i^.(C)l < 1^"^^^''''" 



{dR - 1)9(Z) • 

The rest is simply a calculation of Cl and Or. D 

Now, we take a sequence of graphs as above and let ^ grow to infinity. 

Lemma 16. Let G be a sequence of deterministic graphs, each with a root with an acyclic 
2^ + 1 neighborhood. Suppose d/j — )• oo and C — ^ oo. Fix e > 0, and let z be a sequence with 
|Re(z)| > e and \Q{z)\ > l/dR. Suppose that r^^^ = o(l/d^). Then either 

{{dR - l)-^/'A - .);;^^ _^ - sl{z)\ = 0,{l/dR), 

or 

{{dR - ly^i'A - z);;„^ _^ - sr{z)\ = o,{i/dR), 

depending on whether the root of G has degree dL or dR. We use Oe(-) to indicate that the 
constant in the expression depends on e. 

Proof. Consider the case where the root of G has degree dL- We will proceed in two steps, 
bounding first |(^l(2C) — sl{z)\ and then El{2C,). 
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We define the quantity 



/3 



q2(;{z) - (dR-l) ^q2(:-i{z) 
1 

l-{dR-l) ^ l_^-2C-2 



This is of interest because ipii'^C) = /^^lC^C)- With the assumptions of this lemma, one can 
calculate directly that |/3| = 1 + 0{l/dR); by these assumptions and Lemma [T^ \ipi{2() — 
sl{z)\ = Oe{l/dfi). Since sl{z) is bounded for \z\ > e, this also implies that (PlC^C) = Oe(l)- 
Thus 

1^20 -SL{z)\ = \f^ipL{20-SL{z)\ 

<\^L{2C)-SL{z)\+\i/3-l)ipL{2C)\ 

< 0,{l/dR) + 0,{l/dR) = 0,{l/dR). 

Since tPl{2C) = /3V'l(2C)) by Lemma [HI and our bound on /3, 

\iiLm\ = a'^/\dR-ir'^o,ii/dR). 

Combining this with our bound on El{2C,) from Lemma [T5] gives 

These two bounds prove the lemma. The case when the root of G has degree dR is the 
same. D 

4.4 From deterministic graphs to random graphs 

The main actors of this section will be sequences s, C, and r/. We will choose s and z in 
such a way that s > r[z), and C, will represent the size of an acyclic neighborhood in the 
graph. We will choose r] so that we can control the Stieltjes transform of our graph on the 
set U = {z : 9(2;) > /?}. The following lemma gives us the relation between s and r/: 

Lemma 17. Let r = r{z), fix some s > 1, and let r] = s^'^ — s~^'^. IfQ{z) > r], then r > s. 

Proof. First we prove this when a = 1. Consider the set Eg = {z : \F{z)\ < s}. This set is 
the interior of an ellipse whose foci are —1 and 1 and whose radii are 2(s + 'S~^) and 2{s — s~^) 
(see |MH031 p. 14]). It suffices to show that if Q{z) > ■>], then ^z^ — 1 lies outside of Es- To 
this end, we note that the transformation z 1— )■ ^z"^ — 1 takes the region given by Q{z) > r] to 
the region bounded on the right by the parabola P = {^(t^ — ry^) — 1 + rjti : t G M}. This can 
be checked to touch Eg at — 2('5 + s~^) and otherwise to lie to the left of it. This proves the 
lemma when a = 1. 

To extend this to the case where a > 1, we consider the image of Es under the map 

z I— 7- (a~^'^(22; + a^'^ + a^^'^)) , which is the inverse of 2-v/a(-2^ ~ ck~^ ~ !)• Our argument 
for a = 1 establishes that in this case, the maximum imaginary part of this set is rj. It is 
straightforward to check that for any z, the quantity 9(a^^'^(2z + a"^'^ + a~^'^)) decreases 
as a increases, which establishes the lemma. D 
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For the remainder of this section, let G be a random biregular bipartite graph on m + n 
vertices satisfying ([I])-([3]) as well as the condition cIr = exp (o(l)-v/Togri). Let A be the 
adjacency matrix of G. We define the sequences 



a = mm 



logn 



9(logd,j)' 



7,dR] , 



s = e 



l/a 



logn 
' log dR 



1, 



We now show that sufficiently many vertices of G have tree-like neighborhoods. 
Lemma 18. Let J be the set of vertices in G whose 2(^ -neighborhoods are acyclic. Then 



1 



1^1 > V 



n + m dR 



o(l/n). 



Proof. This is nearly the same as Lemma [5l We may assume C, is an integer by replacing it 
with [CJ . We define 

2C 



N*=Y,HdR-l?'^-'X, 



i=2 



recalling that Xi is the random variable denoting the number of 2i-cycles in G. We have the 
bound n + m — \J\ < N*. Now we apply Proposition H] to calculate 



2C 



E[iVl = Yl 2«(rffi - 1)'^~V* 1 + O 



i=2 



i{i + dR) 



n 



o{d'k) 



and 



2C 



Var[iV*] < 2CY,^e{dR - lf^-^%i 1 + 



i=2 

c 



d%{ia^'-' +a-'dR) 



n 



<2CY,'^i{dR- 1^^(1 + 
By Markov's inequality, 



n'^^Ca^^-^ + dR) 



n 



\J\ > V 



n + m dR 



<P 



< 



N* > 



(1 + a)nr] 
dR 

o{dr'+edr') 



< 0{n-^d-^* + n-3/2^3-) ^ ^^/^), 



D 
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Let Jl and Jr denote the sets of vertices with acychc 2C-neighborhoods in the left and 
right vertex classes, respectively. As in Corollary [6l it is immediate that 



m- \Jl\ ^ r/ 



n + m 



dR 



o(l/n) 



and 



n- \Jr\ ^ _ri_ 
n + m ~ dfi 



o(l/n). 



It is also straightforward to see that this lemma holds when we require the vertices to have 
acyclic 2(^ + 1 neighborhoods rather than just 2^ neighborhoods. 

We can now apply all of these results to the task of bounding the rate of convergence of 
the Stieltjes transform: 

Theorem 19. Fix some e and let U denote the set of complex numbers 

U = {zGC: Re(z) > e, 9(z) > ??}. 
Let Sn{z) denote the Stieltjes transform of {dji — 1)^^^. Then for sufficiently large C^, 



sup \Sn{z) - S{z)\ > Ce/dR 



o(l/n). 



Proof. Let J denote the vertices with acyclic 2C + 1 neighborhoods. We condition on the 
event that (tti— | Ji|)/(m + n) < r]/dR and (n— | J/j|/(m + n) < r]/dR, which by the discussion 
following Lemma [18] holds with probability 1 — o(l/n). We call this event 0. We compute 
Sn{z) for z ^ U, breaking up the vertices into Jl, Jr, and the remaining vertices: 



+ ^Y.i(dR-l)-'/'A-z);l. 
We begin with the third term. Applying the bound ((d/j — \)~^''^A — z) < r/"^, we have 



^—Y.{^dR-l)-'/^A 

1_ rn Z — J ^ 



n + m 



V^J 



m + n —\J\ 1 

< ■ — - < — 

(m + n)ri dR 



(16) 



on the event 0. 

Since every vertex in J has an acyclic 2^ + 1 neighborhood, we will apply Lemma [16] to 
estimate the first two terms. First, we confirm that the conditions of the lemma hold. By 
expanding t; as a power series, we deduce the bound rj > 1/a > 1/dR. We can calculate 



-2C 



exp — 



logn 
a ySlogdR 
9 



1 



<exp(-^log(iR + o(l)) =0(1/4). 
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By Lemma [T71 it follows from 9(z) > i] that r{z) > s. Hence for any sequence z £ U, we 
have r{z)~'^'^ < s~^^ = o(l/d'j^). Thus the conditfons of Lemma [161 hold . and so for all v G Jl, 



and for all v £ Jr, 



{{dn - l)-'/'A - z);l - sdz) = 0,{l/dn), 



{{dR - l)-'/^A -z)^l- sr{z) = 0,{l/dn). 



Combining these estimates with (J16p . 



Sniz) 



siiz) + 



\Jr\ 



■SR{z)+0,{l/dR) 



m + n ^ ' m -\- n 

TTi ft 

= —sl{z) + ^sr{z) + 0,{l/dR) 

m + n 771 + n 

= s{z) + 0,{l/dR) 
on the event fi, which proves the theorem. 

We now restate and prove our local convergence law. 
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Theorem 3. Fix e > 0. Let ^„ be the BSD of (dR — 1)^^/^A, and let /i be the limiting BSD 
defined in Corollary \^ There exists a constant C^ such that for all sufficiently large n and 
5 > 0, for any interval /CM avoiding [— e,e] and with length \I\ > max (irj, r] / (—6 log 6)^ , it 
holds that 



\fIn{I)-Kl)\ <SC,\I\ 



with probability 1 — o(l/n). 



Proof. This theorem follows from the arguments of [TVIH Lemma 64], which we will sketch. 
Define 



ny) = - [ 



rp' + {y — xY 



dx. 



This function F{y) approximates the indicator function on the interval /, and the following 
statements hold: 



j F{y) dfi{y) = fi{I) + O, L log ^"j , 



/ 



F{y) d/i„(y) = fin{I) + Oe ( 7?log 



The proofs of these statements in [TVIH Lemma 64] have /i as the semicircle law, but they 
apply just as well to our limiting measure fi; the only thing necessary for the proof to go 
through is that /x has a bounded density outside of the interval [— e, e]. On the event Q of the 
previous theorem. 



Fiy)dfiiy)- J F{y)df,niy) 



< 



1 

vr 

C.|/| 

irdR 



{x + r]i)) — Q{sn{x + r]i))) dx 
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As observed in [TVllj . it follows from the condition |/| > r]/{—6log6) that r^log '-^ = 0{6\I\). 
Since d_R — )• oo, for n sufficiently large, 

|/i(/)-/i„(/)|<C7e5|/| 

for some constant Ce (not necessarily the same one as before) on the event i}. D 

Appendix 

Our goal here is to prove Proposition [H We mention that it is possible to prove much more 
than this. The main result of [MWW04J is that the distribution of short cycles in a random 
regular graph is approximately Poisson, and this result holds for biregular bipartite graphs 
as well, with suitable modifications of the proofs. 

We will use a theorem from [McK81| that gives us the probability that G contains some 
subgraph L C Km,n- For any v € Km,n, let gv and l^ denote the the degree of v considered 
as a vertex of G and of L, respectively. Let /max be the largest value of k. Consider L to be 
a collection of edges, so that \L\ is the number of edges of L. The notation [x]a denotes the 
falling factorial, x{x — 1) • • • (x — a + 1). 

Proposition 20. Let L C Km,n- 

(a) If \L\ + 2dR{dR + /max - 2) < ndn - I, then 

P[L C G] < M^^ — 

(b) IfndR - 2dR{dR + /max - 1) - 1 - |i| > d/j/max, then 
W[9i\h 



P[L C G] > 



[ndR-l]\L 



\L\ 



1 (^fl/max ) / ( 1 _|_ R 



\u 



ndR - \L\ - 2dR{dR + Invax - 1) - 1 J \ ndfl - 2dR(di? + /max - 2) - 1 - |L|(e - l)/ey _ 

Proof. This is an application of |McK8H Theorem 3.5], with the set H from that theorem 
equal to 0. D 

We are most interested in when L is a cycle with 2r edges, in which case the above 
theorem reduces to the following: 

Corollary 21. Let L C Km,n be a cycle with 2r edges whose presence in G we wish to test. 

(a) If 2r + 2(i| < ndR - 1, then 

^ - ^ - [ndR - 24 - 1]2. 

(b) IfndR - 2dR{dR + 1) - 1 - 2r > 2dR, then 
dlidL - lyd'^idR - 1) 



\r 



P[L C G] > 



[ndR - 1] 



2r 



2di? \ l(. , 4 



1 + 



ndR-2dR{dR + l)-l-2r) I \ ndR - 2d%- 1 - 2r{e - l)/e 
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2r 



Proof of Proposition \^ Let L C Km,n be a cycle with 2r edges. We start by showing that 



n^^ V \ n ' ndji 



P[LCG]= '-^ -^^7^ -^- (l + 0(— + r:^))- (17) 



Since d^j = o{n), all of the conditions for Proposition [20] and Corollary [21] apply. By 
Corollary [2iM 



P[L C G] < 



(dL - l)''(di; - lya- '-d^' 
{ndR - 2d\ - 2r) 

{dL - IfidR - Ifa-' ( ndR 



R 

2r 

2r 



n2^ V^^rfi?^ - 2dR - 2r 



(di - lYidR - ly-a-'- (^ , 24 + 2r 



i? 

2r 



1 + 



n^*^ \ \ n ndij ' ' 

The last line follows from the fact that if a; > —1, 

(1 + xf < e^^ = 1 + OM 

as rx —7- 0. 

From Corollar\ t2]lb] 



{ndRf 
1- .. .. ''^. . . 1/(1 + 



ndij - 2dij(dij + 1) - 1 - 2ry ' V ridR -2dj.-l- 2r(e - 1)/ 



2r 



(dL - ir(dR - Ifa-'- / l + 0(l/n) 





„2r 




{dL- 


- l)^(di? - 


- l)^a-^ 




^2r 




{dL- 


- l)'(di? - 


- l)'"a-'" 



2r 



1 + 0{dR/n) 



,i + o 



n^' \ \ n 



n 
rdR 



These two inequalities prove p!7|) . 

The number of cycles of length 2r in Kjn,n is [r7T-]j.[n]r/2r. Using [n]r = n^{l + 0{r'^/n)), 
which can be proven by showing inductively that [n]r > n^(l — r^/2n), we calculate the 
expected number of such cycles: 



2r n^^ \ \ n ndR 

(na)'n'-(l + 0{ryn)) {dL - lYidR - If a-' (^^^ /rdR ^ r^ 



2r n^'^ \ \ n ndR 
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Let C be the set of 2r-cycles in Km,n- We will calculate Var[Xr] using the equation 
We break up C x C to help calculate this sum. 



Ci = {(C7i,C2)eCxC 

C2 = {(Ci,C2)GCxC 
C3 = {(Ci,C2)GCxC 



CinC2 = 0}, 
CinC2/0, butCi^C2}, 
Ci = C2}. 



We are considering cycles as collections of edges, so pairs of cycles that share vertices but 
not edges belong in Ci rather than €2- 

For (Ci,C2) G Ci, Proposition I2Q61 and a calculation identical to the one in (llSh show 
that 



P[Ci U C2 c G] < (^^ - ^^'^(^fr ^^'^""'^ ( 1 + O ( — + 



n' 



4r 



n ndij 



Bounding |Ci| by \C\ = ([m]r.[n]r/2r)^ and repeating the calculation in (fT9]) gives 

r(r + (Jr) 
n 



E P[CiUC2CG] <;u2(^l + 

(Ci,C2)6Ci 



(20) 



For a lower bound on this sum, we note that by Proposition I2db[ for any 2r-cycles Ci and 
C2 that share no vertices, 



P[Ci U C2 c G] > 



di^idL - irdj^idR - ly- f 1 + 0(l/n 



(nd/j) 



4r 



1 + 0(dfl/n) 



4r 



(di - l)''^(d« - l)^^a-^'- ^ + 0^'"'^^ 



n' 



4r 



n 



Summing this over the ([n]j.[?n,]j./2r)([n — r]r[m — r]r/2r) such pairs of 2r-cycles provides the 
bound 



E P[<^i U G2 c G] > 
(Ci,C2)eCi 



[n]r[m]r[n - r]r[m - r]ra ^'' 2 (. . q f^dn 



n 



4r 



n^m''(n — r)'~(?n, — r)'~(l + 0(r^/n))a ^^ 2 f /^ ( ^^R 



n' 



4r 



n 



n) \ an) I ^ \ \ n I I \ \ n 



Mr 1 + O 



r{r + di?) 



n 



The sum over Cs is 



E P[CiUG2CG]=E[Xr] = /i.(l + 0(^^^^^^^ 



(21) 



(Ci,C2)eC3 
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To estimate the sum over C2, we bound the number of isomorphism types of a graph 
H = CiUC2 for (Ci,C2) G C2. Let H' be the graph {V{Ci) D V{C2),E{Ci) D E{C2)). Say 
that H' has p components and j edges. As H' is a forest, it has p + j vertices, so H has 
ir — p — j vertices. We also note that H has 4r — j edges. 

Let Ci = aia2 • • • a2r and C2 = 61^2 • • • b2r, with ai = bi. Let Ai, . . . ,Ap be the compo- 
nents of H' ordered as the appear in Ci. We can encode the isomorphism type of H in the 
following four sequences: 

• Si is the number of vertices in Ai 

• ti is the smallest j such that aj G Ai 

• Ui is the smallest j such that bj G Ai 

• fj specifies whether Ai is oriented the same way in Ci as in C2 (if ^j is a single vertex, 
consider it to be oriented the same) 



For example, the following diagram is the union of two cycles of length 
black and the second colored gray. 



the first colored 




The intersection graph H' has three components, Ai = ai = bi, A2 = 0304 = 6263, and 
A3 = QQa-ras = bQbjbs- The four sequences for this graph are 

s: 1,2,3 
t : 1,3,6 
u: 1,2,6 
V : yes, yes, no 

To illustrate that the isomorphism class of H is encoded in these sequences. We will demon- 
strate how to recover it in this example. Start by drawing a cycle and labeling its vertices 
oi, . . . ,08. From s and t, we can deduce that Ai = oi, A2 = 0304, and A3 = aQUjas- From 
u and V, we deduce that 61 = oi, 62 = 03, 63 = 04, be = as, by = 07, and bs = aQ. Since 64 
and 65 are unaccounted for, we conclude that they are not contained in oi • • • as. Once we 
add edges to connect 61 to 62, &2 to 63, and so on, we have recreated H up to isomorphism. 
Now, we consider the number of possible isomorphism classes of some H = CiL)C2- The 
sequence s is a composition of p + j into exactly p parts, so there are (^^^7 ) possibilities. 
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We know that ti = 1, and we know that t2, ■ ■ ■ ,tp are ordered and are distinct, so there are 
at most ( ''"-j^) choices for t. We know that ui = 1 and that U2,- ■ ■ ,Up are distinct but not 

necessarily in any order, so there are at most (p~i)(p ~ 1)' choices for u. For v, there are 
2^^^ choices. For a fixed choice of of p and j, the number of possible isomorphism classes of 
H is hence bounded by 

By replacing p + j — 1 and 2r — 1 by 2r, we bound this quantity by 

Suppose an isomorphism type has a vertices from the left vertex class and b from the 
right vertex class (with a + h = Ar — p — j). Then this isomorphism type can be realized in 
at most [n]a[m]b + [n]b[m]a < 20^^'"^'^ n^^~^^~^ ways. By Proposition [2Qfal the probability of 
any one of these realizations being a subgraph of G is bounded by 

[ndR-U\-l]ir-j ~ a(4'^-i)/2n4^-i ^^'' 

All together, we have the bound 






Combining this with (j20p and (j2ip shows 

- ^[Xrf + O 



n / / \ \ n 



n 






^r + fJ'rO 



r{a(^-'^y^r + dR) 
n 



Hence 



Var[X^] = /i J 1 + ^^O ^ ^ 



n 



M. 11 +01''"'^"""^°"''"' 11. 



n 



27 



References 

[BL13] Shimon Brooks and Elon Lindenstrauss. Non-localization of eigenfunctions on 

large regular graphs. Israel Journal of Mathematics, 193(1):1-14, 2013. 

[BSIO] Zhidong Bai and Jack W. Silverstein. Spectral analysis of large dimensional ran- 

dom matrices. Springer Series in Statistics. Springer, New York, second edition, 
2010. 

[DJPP12] loana Dumitriu, Tobias Johnson, Soumik Pal, and Elliot Paquette. Functional 
limit theorems for random regular graphs. Probability Theory and Related Fields, 
pages 1-55, 2012. Published onhne, 25 August 2012. 

[DP12] loana Dumitriu and Soumik Pal. Sparse regular random graphs: Spectral density 
and eigenvectors. Ann. Probab., 40(5):2197-2235, 2012. 

[Dum03] loana Dumitriu. Eigenvalue Statistics for Beta-Ensembles. PhD thesis, Mas- 
sachussetts Institute of Technology, 2003. 

[EPR+10] Laszlo Erdos, Sandrine Peche, Jose A. Ramirez, Benjamin Schlein, and Horng- 
Tzer Yau. Bulk universality for Wigner matrices. Comm. Pure Appl. Math., 
63(7):895-925, 2010. 

[ERS^IO] Laszlo Erdos, Jose Ramirez, Benjamin Schlein, Terence Tao, Van Vu, and Horng- 
Tzer Yau. Bulk universality for Wigner Hermitian matrices with subexponential 
decay. Math. Res. Lett, 17(4):667-674, 2010. 

[ERSYIO] Laszlo Erdos, Jose A. Ramirez, Benjamin Schlein, and Horng-Tzer Yau. Uni- 
versality of sine-kernel for Wigner matrices with a small Gaussian perturbation. 
Electron. J. Probab., 15mo. 18, 526-603, 2010. 

[ESY09a] Laszlo Erdos, Benjamin Schlein, and Horng-Tzer Yau. Local semicircle law 
and complete delocalization for Wigner random matrices. Comm. Math. Phys., 
287(2):641-655, 2009. 

[ESY09b] Laszlo Erdos, Benjamin Schlein, and Horng-Tzer Yau. Semicircle law on short 
scales and delocalization of eigenvectors for Wigner random matrices. Ann. 
Probab., 37(3):815-852, 2009. 

[GM88] C. D. Godsil and B. Mohar. Walk generating functions and spectral measures 
of infinite graphs. In Proceedings of the Victoria Conference on Combinatorial 
Matrix Analysis (Victoria, BC, 1987), volume 107, pages 191-206, 1988. 

[Johl2] Tobias Johnson. Exchangeable pairs, switchings, and random regular graphs. 
Preprint. Available at arXiv:1112.0704, 2012. 

[JP12] Tobias Johnson and Soumik Pal. Cycles and eigenvalues of sequentially growing 

random regular graphs. Preprint. Available at arXiv:1203.1113, 2012. 

[LW08] Lisa Lorentzen and Haakon Waadeland. Continued fractions. Vol. 1, volume 1 
of Atlantis Studies in Mathematics for Engineering and Science. Atlantis Press, 
Paris, second edition, 2008. 

28 



[McK81] Brendan D. McKay. Subgraphs of random graphs with specified degrees. In 
Proceedings of the Twelfth Southeastern Conference on Combinatorics, Craph 
Theory and Computing, Vol. II (Baton Rouge, La., 1981), volume 33, pages 213- 
223, 1981. 

[MH03] J. C. Mason and D. C. Handscomb. Chebyshev polynomials. Chapman & 
Hah/CRC, Boca Raton, FL, 2003. 

[MS03] Hirobumi Mizuno and Iwao Sato. The semicircle law for semiregular bipartite 
graphs. J. Combin. Theory Ser. A, 101 (2): 174-190, 2003. 

[MWW04] Brendan D. McKay, Nicholas C. Wormald, and Beata Wysocka. Short cycles in 
random regular graphs. Electron. J. Combin., ll(l):Research Paper 66, 12 pp. 
(electronic), 2004. 

[Sta99] Richard P. Stanley. Enumerative combinatorics. Vol. 2, volume 62 of Cambridge 
Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 1999. 

[TVIO] Terence Tao and Van Vu. Random matrices: universality of local eigenvalue 
statistics up to the edge. Comm. Math. Phys., 298(2) :549-572, 2010. 

[TVll] Terence Tao and Van Vu. Random matrices: universality of local eigenvalue 
statistics. Acta Math., 206(l):127-204, 2011. 

[TVW13] Linh V. Tran, Van H. Vu, and Ke Wang. Sparse random graphs: Eigenvalues and 
eigenvectors. Random Structures Algorithms, 42(1):110-134, 2013. 



29 



